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Abstract 

We study the role of the Jacobian arising from a constraint enforcing the 
nonlinear relation: g = pV, where p, g and V are the mass density, the 
momentum density and the local velocity field, respectively, in the field theo- 
retic formulation of the nonlinear fluctuating hydrodynamics of simple fluids. 
By investigating the Jacobian directly and by developing a field theoretic 
formulation without the constraint, we find that no changes in dynamics re- 
sult as compared to the previous formulation developed by Das and Mazenko 
(DM). In particular, the cutoff mechanism discovered by DM is shown to be 
a consequence of the 1/p nonlinearity in the problem not of the constraint. 
The consequences of this result for the static properties of the system is also 
discussed. 
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I. INTRODUCTION 



The development of the appropriate field theoretical treatment for nonlinear fluctuating 
hydrodynamics (NFH) of simple fluids is a more subtle enterprise than one might first 
imagine. If one develops a theory for simple fluids including the complete set of conserved 
fields, the mass density p(x), the momentum density g(x), and the energy density e(x), 
then one finds M that one must include multiplicative noise in order to gain consistency 
with thermodynamics. This surprising result is associated with the connection between the 
fluctuating energy, entropy and temperature. In the case where the energy is not included in 
the set of slow variables, one still finds technical problems in developing the associated field 
theory. These problems are associated with the form of the kinetic energy density g 2 /2p 
resulting from Galilean invariance. This 1/p factor can be identified with the nonlinear 
relationship between the momentum density and the local velocity field: V(x) = g(x)/p(x) 
and, at first sight, seems to complicate the problem considerably. 

In Ref. (DM), the NFH of compressible fluids was studied as a model for the glass 
transition. A field theoretic formulation of the problem was developed by generalizing the 
standard Martin-Siggia-Rose (MSR) method [|J to include the nonlinear constraint between 
g(x) and V(x). In the functional integral formulation described in Ref. [[J, this constraint 
was enforced by introducing an auxiliary velocity field V(x) and inserting a delta-functional 
constraint, 

J PV(x)5(g-pV). (1) 

Enforcing the relation g = pV eliminates the 1/p nonlinearity in the kinetic energy and 
yields a polynomial action in p(x),g(x) and V(x), which in turn allowed one to carry out 
a perturbation theory expansion with standard renormalization schemes. The effects due 
to nonlinearities on various physical quantities were then calculated at the one-loop order 
and, as a result of a systematic perturbative expansion, it was discovered that there is a 
nonhydrodynamic correction that cuts off the sharp nature of the idealized glass transition 
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H|7|]. This result has important consequences that dense fluid systems remain ergodic for all 
values of the density and temperature, although the density feedback mechnism does drive 
the viscosity to large values. 

Recently Schmitz, Dufty and De || suggested that a more faithful representation of the 
original Langevin equations requires that the constraint should be enforced by using, instead 
of Eq. (|), the identity 

1 = / DV(x) S(- - V) (2) 
J p 

= \\p\\ J PV(x)5(g-pV), (3) 

where \\p\\ is the Jacobian resulting from the change of variable in the delta-functional. It 
was argued in Ref. || that the cutoff mechanism found by DM may be an artifact resulting 
from their use of the constraint, Eq. (|l|) rather than Eq. (|2]). Since these constraints differ 
only by the Jacobian factor this is equivalent to saying that the cutoff mechanism is 
eliminated if one includes the Jacobian factor into the development. We show here that this 
is not the case. 

We can evaluate the Jacobian \\p\\ directly by expressing it in terms of an integral over 
Grassmann fields ||. Inserting this result into the DM action, we will be able to show that 
it does not play any significant role on the dynamics. We then develop a similiar functional 
integral formulation to the one in Ref. ||, but without introducing any constraint. The 
1/p nonlinearities are expanded as a power series in 5p = p — po, and by comparing our 
perturbation theory results with those of Ref. 0, we will be able to see explicitly the role of 
the Jacobian. We find that the two formulations are equivalent at the one-loop order of the 
perturbation theory and the Jacobian has no dynamical effects at this order. This seems 
perplexing since it is clear that the Jacobian does influence the static equilibrium behavior 
generated by a general effective Hamiltonian. The resolution to this apprent contradiction 
is that the dynamics of a simple fluid, generated by a Langevin equation, is insensitive to 
changes in the chemical potential. It appears, at least to lowest order, that the effects of the 
Jacobian in question can be absorbed in the chemical potential and does not directly affect 
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the statics generated by the Langevin equation. 

In Sec. m we give a brief review of the formulation in Ref. || just to collect some results 
needed for our discussion. In Sec. [TTJ, the Jacobian is expressed as a Grassmann integral and 
the usefulness of this formulation is discussed. In Sec. [TV], the theory is formulated without 
introducing the velocity field and the one-loop equivalence between the two formulations are 
studied in detail. The static theory is described in Sec. |V|. 



II. BRIEF REVIEW OF DM 

Our starting point is the set of generalized Langevin equations for compressible fluids 
for the set of slow variables {p(x),g(x)}. Following standard procedure ||10|| , we obtain the 
continuity equation for the conservation of mass 

(4) 

and the generalized Navier-Stokes equation for the conservation of momentum, 

f = -pv^ - £ v^-f) - E LlJ ( 9 j) + e, (5) 

In Eq. (|5|), F u [Sp] is the potential energy part of the effective Hamiltonian given by 

F = F K + F u , (6) 

where the kinetic energy Fk is 

f> = /Af!. (7) 

In Ref |H, the simple choice of F u , 

F u [5p] = J d d ^{8pf (8) 

was used, where A^ 1 = po/c^ is the flat static structure factor, p is the average mass 
density and Co is the bare sound speed. This simple quadratic part of F u then gives the 
pressure nonlinearity in the Langevin equation which is responsible for the density feedback 



mechanism of the glass transition ||. In general, F u can be any local functional of Sp and 
the spatial derivatives of 5p. By including the derivatives of dp, we can probe the effect of 
the spatial correlations in the system on the dynamics. The dissipative matrix in Eq. (|5]) is 
given by 

L y (x) = -Tjo^ViVj + <%V 2 ) - CoViV,-, (9) 

where rjo is the bare shear viscosity and Co the bare bulk viscosity. For later use we define 
the bare longitudinal viscosity r = Co + f^o- The noise 0j is Gaussian with variance 

(6i(x, t)Qj(x', t')) = 2k B TL ij (x)5(x - x')5(t - t'). (10) 

These Langevin equations can be put into a field theoretical form following the standard 
MSR procedure . It essentially amounts to introducing a hatted variable ip for each field 
-0 = p, g to enforce the equation of motion and integrating over the Gaussian noise to yield 
a quadratic action in ■0. The generating functional without source terms is given by 

Vi/>T>$e- 8 W'®, (11) 

where the action S[ip,ip], with the notation 1 = (xi,£i), is given by 

dp 

giP L,ij{i)gj-t-ip[ 

+ pV{ ^ + S v - ( ^ i) + E )] }• (12) 
i at i °p j p j p 

In Ref. [H, the 1/p nonlinearities are removed by using Eq. ([!]) and the representation of a 
delta-functional, 

S(S ~ PV) = I X?V exp[z | dl V(l) ■ (g(l) - p(l)V(l))], (13) 

so that the generating functional in Ref. || can be expressed as the functional integral over 
\l/ and where \1> = p, g i} V^. Denoting hereafter the results from Ref. || by the supersript 
I, we have 

Z {i) = J Mp^-smif,^ ( 14 ) 



S[^] = J dl{Y / g i r 1 L ij (l)g j + ip[^- + V 1 - S } 

in 1 



where S^[\l/, is given by 

ij 1 
* 1 " 3 3 

+ *V-(g-pV)}. (15) 

The perturbative expansion for a polynomial action of this form is standard. The non- 
linear corrections are given in the form of the self-energies modifying the inverse linear 
propagator through Dyson's equation, 

where G a p = {^f a ^p). Detailed calculations of the self-energies were carried out in Ref. [[|. 
Here we list the results for the density response and the correlation functions. In the Fourier- 
transformed space, 

G ™ (q ' u > ~ DW(q, W ) ' ( ?) 

G "" (q ' u) ~ ^>mp ' (18) 

where 



p(q, u) = po - z"S^(q, w), (19) 
L(q,^) = g 2 r + z4?(q,^), (20) 
L(q,u,) = g 2 r - (2/3- 1 )~ 1 {Sg(q,a;) + (^) 2 E& (q, (21) 



and 



L>« (q, w ) = p(q, cu){^ 2 - g 2 c 2 - E^q, w)} + zX(q, cu){cu + ig£$(q, w)}. (22) 

Here all the self-energies are their longitudinal parts. The density feedback mechanism is 
realized by calculating the one-loop self-energies contributing to the dynamic viscosities, 
L(q, uj) and L(q, u), which are |§ quadratic in the density correlation function [fj. Without 
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the nonhydrodynamic correction due to the self-energy (q, u), the response function is 
the same as the one that gives the density feedback mechanism. Thus S^(q, u) provides 
the DM cutoff. As mentioned in Sec. | and as indicated by the subscripts, the cutoff is 
related to the 1/p nonlinearities and the constraint between g and V. 



III. EVALUATION OF THE JACOBIAN 

The question arises: what changes result if one uses the constraint Eq. @ rather than 
Eq. (jl]) to introduce the velocity field? It should be clear from Eq. (||) that the only changes 
in the DM action comes from the Jacobian factor ||p||. This quantity can be represented as 

||p|| = J V V Vf] exp[/ dl 5>(l)p(lMl)], (23) 

i 

where fji and rji are Grassmann fields. We must, as indicated below, be careful in using 
any unregularized representation of the Jacobian like Eq. ( p3[) in the functional integral 
formalism, since it might not preserve causality. If we include this in the DM formulation, 
then the appropriate generating functional is given by 

Z = J Vr]Vf)V^Vi> e- §[ *'*>™}^ (24) 
where the action is given now by 

§[%%r,,fi] = S^[%9] + J dl (25) 

i 

The only linear propagator involving fj or rj is 

G lvM^) = - G lnMi") = J Q ( 26 ) 

The new nonlinear term arising from the Grassmann fields is of the form, SpJ^iViVi- This 
gives, at one-loop order, only two new self-energy diagrams (Fig. 1). The first appears to give 
a contribution to S pp , which clearly violates causality. We note that the self-energies between 
two unhatted variables are indeed equal to zero JTT| in the original DM formulation in 
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accordance with the causality. Thus Eq. (|23|) has to be regularized. We note that if we use, 
instead of Eq. (^) , 

/ dl £^(l){e^+p(l)Ml), (27) 

with a time derivative with a very small coeffcient e, then the first diagram of Fig. 1 vanishes, 
since it is in the form of a time integral of the product of a retarded and an advanced 
propagator. Furthermore, as will be clear later in the section, the regularization given as 
Eq. (|27|) guarantees the causality to all orders of perturbation theory, i.e. E^ remains to 
vanish to all orders. 

The second diagram in Fig. 1, however, does give a finite contribution. In the limit 
e — > 0, we have 

(q> <*0 = -S»ww (<1> u ) = 5 a(0) = <^E, (28) 

c o 

where 5a (0) = J A -A^s with the large momentum cutoff A. This self-energy gives a correction 
to the correlation function between rj and fj; 

<Wq,o,) = ^-. (29) 
Po — ^ 

We note that there are no such self-energies at one-loop order that link the Grassmann fields 
to the p, g or V. This fact, together with the vanishing of the first diagram in Fig. 1, shows 
that the inverse propagator matrix in Eq. (|T6|) is in a block-diagonal form with the block 
being identical to those given in DM formulation. Inverting the inverse propagator, we find 
that there are no changes, at one-loop order, in the response and correlation functions for 
p, g and V due to the Jacobian. 

This result, in fact, can be generalized to all orders of perturbation theory by deriving the 
Ward identity from the invariance of the action, Eq. ( p5|) under rji — > e la r]i and f/« — > e~ la f]i 
for any constant a. Infinitesimally this symmetry transformation reads 

Srji = {ia)r]i, 5f]i = (-ia)fji, S^ Q = S^ a = 0. (30) 
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Introducing the sources J a , J a for ^/ a , ty a and the Grassmann sources ^ for r]i, 77^ respec- 
tively into the generating functional, Eq. fl2~4]) , we have 



Z[J,J,£,f] = | VrjVfjV^vif exp[ 77, 77], 

+ [ dl[^{JA + iA} + E{^ + ^}]]. (31) 

Since the action 5 and the integration measure in Eq. ( |5T] ) are invariant under the trasfor- 
mation, Eq. fl3U|), the variation of the source terms must vanish if we change the integration 
variable by Eq. (RDT) . This yields the following Ward identity: 



/1 Ete(i)^ ir 6(i)^ iy }^.-'.f.a=o (32) 

Taking a derivative of Eq. (|32[) with respect to £j(3) and then with respect to J a {2) and 
setting all the sources equal to zero, we have 



= -Z[0] 0*^.(2, 3). (33) 

J=J=£=£=0 



Similarly, we can easily derive 

= = Gy- = 0. (34) 

It follows that, as in the one- loop case, the inverse propagator is block-diagonal. Let us 
now investigate whether the nonlinearity, 5p Y^i ViVi contributes to the self-energies in the 
W block £^ and S^,^,). To do this, we first note that Eqs. ([33]) and ( |34"D 

can be understood as a direct consequence of the charge conservation for rj, f] imposed by 
the symmetry, Eq. (0). We can imagine a charge flowing through the rjfj lines in the 
perturbative diagrams, say from rj to fj. The vanishing of the correlation functions with only 
one external rj or fj, Eqs. ( |33| ) and (|34]), then directly follows from the charge conservation. 
The internal rjfj lines must then form a complete circle following the charge flow. The most 
general diagram involving rj, fj that can contribute to the self-energies in the tyi£f block is 
drawn in Fig. 2. At one-loop order, it is just the first diagram in Fig. 1, which vanishes JTT 



because of the regularization, Eq. (|27|) . We have, in fact, the same situation for the general 



diagram in Fig. 2. This diagram is expressed in terms of the integral over I loop momenta 
and frequencies, where I is the number of the loops in the diagram. Now we inverse Fourier 
transform the zeroth order r)fj propagators to get the time integral over t\,---,t n , where n 
is the number of G°=(ij)'s in the diagram. Then we can integrate over one frequency, which 
yields the delta function S(J2iti), therefore we have the following factor inside the integral 
over / momenta and I — 1 frequencies: 

„ n n 

j Rdt i 5{Y J U)G% i {t 1 )---G% i {t n ) = U, (35) 



i=l i=l 

which vanishes due to the regularization, Eq. (|27|), since all the G^(tj)'s are advanced 
propagators. This shows that the nonlinearity involving the Grassmann fields does not affect 
the self-energies in the block, especially = even after including the regularized 
Jacobian Eq. fl27|) . We can therefore conclude that the Jacobian, included in the DM action 
as an integral over the Grassmann fields, as in Eq. (f25|), has no effect on the correlation and 
response functions for p, g and V at all orders of perturbation theory. 

IV. FORMULATION WITHOUT V FIELDS 

Schmitz et al. || have suggested that the cutoff mechanism found by DM may be some- 
how introduced artificially through the implemetation of the constraint condition. By look- 
ing at the original theory without the constraint within perturbation theory, we show that 
this worry is without foundation. Starting from the original action Eq. (|T2"D, we expand 1/p 
as a power series in 5p. Keeping nonlinear terms that are relevant to the one-loop order, we 
have in the action, Eq. (|T2"D, 

E V 3 -(^).EV J (^)-EV,(^) + -, (36) 

E Mf ) - E M*) - E + EM#) + - (a^) 

j P j P° j Po j Po 

Let us denote the results of this formulation by the superscript II. Only three kinds of self- 
energies are generated by these nonlinearities: £^P(q, uj), E^(q, u) and E^(q, lo). Using 
these self-energies we can represent various correlation and response functions. For example, 
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^ faM) = ggg^z5^0} , (39) 



|fl<">(q,w)| 

where 



Z^)(q,c) = W 2 - g 2 c 2 - Ef/)(q,a;) + £><?u> - ^(q.w). (40) 

Po 



We note that the number of self-energies is reduced from seven to three compared to the 
previous case, since only two variables, p and g, are considered here. But the number of 
diagrams we have to consider for each self-energy is increased according to the appearance of 



the new nonlinearities in Eqs. ( pq) and (37). At one-loop order, by comparing the nonlinear 



vertex structures of the two actions, Eq. (IT4 ) with the expansion Eqs. (pq) , (p7|) and Eq. (|15|), 
and by using the relations among the linear propagators: 

G^(q,o;) = ^G^(q,u;), (41) 



G^(q^) = ^G^(q,^), (42) 



with T being replaced by r/ for transverse propagators, we find that the new diagrams for 
each self-energy of (II) are in one to one correspondence to the diagrams of the corresponding 
self-energies of (I) that are not present in formulation (II). For example, in Fig. 3, the first 
five diagrams contributing to reproduce the corresponding diagrams of S^. The 

remaining diagrams contain the ggp vertex from Eq. ([37]) with g on external legs, which 
give us a factor of for each diagram. Other than this factor, we find that these are the 



a. 2 

PO 



diagrams of the self-energy £^ . Thus, at one-loop order, we have 



= 4 j pW) + 4>^)- ( 44 ) 



Similarly, as seen from Figs. 4 and 5, we find that 



4 7 /W) = -£#(q,w) + (q, W ), (45) 
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S«(q, W ) + (^) J E^(q, W ), (46) 



since £^(q, u;) = — E^(q, uj). Using these relations we can express the response and 
correlation functions in formulation (II), Eqs. (|38|) and (^), in terms of the self-energies in 
formulation (I). We first note that, since we are dealing with one-loop corrections, we can 
neglect quadratic and higher order terms in E's. Inserting Eqs. ([44]) and ((45|) into Eq. (p0|). 
we have 

D( n )(q, W )=p- 1 (q lW )D( : )(q lW ). (47) 
From Eqs. (|38|), fl39|) and using Eqs. (p4|)-(f46|) and (|47|), we have, at one-loop order, 

G%\w) = Gg(q,w), G^\^u) = Gg?(q,uO, (48) 
which shows that the two formulations, (I) and (II) are equivalent to first order. In particular, 



the DM cutoff is recovered as seen from the second term in the right hand side of Eq. ([44]) 
and Eq. ([47|) . Therefore the cutoff is generated not by a particular form of a constraint 
condition, but by the intrinsic 1/p nonlinearity in the problem. 

The renormalizations of various hydrodynamic quantities for the two formulations are 
also equivalent. For example, the renormalized sound speed, for both (I) and (II), is given 
by 

c 2 = c 2 + lim -£^(q,u;). (49) 

q,(jJ^V q 

But, because of Eq. ([P]), limi£^ = limiS^. This quantity was evaluated in Ref. J3J to 
show that 

q,u^0 q yr 

i.e., there is no one- loop correction to the hydrodynamic sound speed. 
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V. STATIC THEORY 



We have investigated in Sec. and Sec. |V| the role of the Jacobian in the dynamics 
governed by the Langevin equations, and found, to lowest order, that there are no changes 
in the dynamics due to the Jacobian. Let us now look at the static equilibrium behavior. 
The static equal time averages of the fields ipi are calculated with respect to the effective 
Hamiltonian, Eq. (|6|): 

(ViVi> = / Ttye-^ify/Z, (51) 

where Z = J Dip e~ 1317 '. Nonlinear corrections come from the 1/p factor in the kinetic energy, 
Eq. (|7|). We note that the potential energy part of the effective Hamiltonian F u appears in 
the Langevin equation only in the form of pV*(5-F u /5p). Therefore a chemical potential-like 
term in F u , which is linear in p, does not affect the Langevin equation. But, in the statics, 
such a term controls the renormalizations of the average (p) of the mass density and the 
sound speed c. To be specific, let us take the simple example 

F u = J d*x[^f-w], (52) 

which gives the same Langevin equation as Eq. (|8|). The one-loop (i.e. 0(/5 -1 )) corrections 
to (p) can be calculated by expanding p = (p) + Ap and integrating out the g field. As a 
result, we have the following one- loop effective action in Ap: 

+{A(p)-^-r 1 ^}(Ap). (53) 
Then (p) is determined by setting the linear term in Ap to zero. Thus we have 

(p) = p + rVi + o(r 2 ), (54) 
pl ~A + ^Ap-;' (55) 

where p = Ap + (3~ l p\ + 0(/3~ 2 ). The sound speed is given by the usual thermodynamic 
result, 
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c 2 =/3- 1 (p)(A P Ap)- 1 . (56) 

The derivation of Eq. ([56]) from the Langevin equation follows using the Fokker-Planck 
description of the problem ||12|| . The one- loop correction to the static density- density cor- 
relation function can easily be read off from the quadratic term of the effective action, 
Eq. (0): 

^(ApAp) = i-/?- 1 ^| + 0(r 2 ). (57) 

Therefore, 

c 2 = c 2 + /r 1 ( Ml + — ) + o(/r 2 )- (58) 

Po 

This is consistent with the dynamics result, Eq. (|50|) , if we choose /ii = — (15^(0) /p . Thus, 
in some sense, the system selects the corresponding static limit as it evolves through the 
Langevin equations. The above discussion indicates that any effect the Jacobian might have 
on the static renormalization can be absorbed in the chemical potential fi chosen by the 
Langevin equation. 



VI. CONCLUSION 

We have shown that the Jacobian, Eq. (|3|) arising from the constraint enforcing the 
nonlinear relation between g and V does not change the dynamics of Das and Mazenko. In 
particular, the functional integral formulations with and without a constraint are shown to 
be equivalent at least to lowest order. Therefore the cutoff mechanism discovered in Ref. || 
is a genuine effect of the 1/p nonlinearity in the problem and not an artifact of the particular 
form of the constraint condition. 
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FIGURES 

FIG. 1. The one- loop diagrams generated by the nonlinearity due to the Jacobian. The Jaco- 
bian is regularized such that the first diagram vanishes. 

FIG. 2. The most general diagram generated by the nonlinearity due to the Jacobian without 
external Grassmann fields. The arrow indicates the charge flow. 

FIG. 3. The detailed correspondence between the one-loop diagrams contributing to E^i and 

FIG. 4. The detailed correspondence between the one-loop diagrams contributing to the 
self-energies in Eq. (|45|). 

FIG. 5. The detailed correspondence between the one-loop diagrams contributing to the 
self-energies in Eq. (|46|). 
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